Thermal noise suppression: How much does it cost? 
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qq , In order to stabilize the behavior of noisy systems, confining it around a desirable state, an 

f"^ ■ effort is required to suppress the intrinsic noise. This noise suppression task entails a cost. For 

the important case of thermal noise in an overdamped system, we show that the minimum cost 
£SJ , is achieved when the system control parameters are held constant: any additional deterministic or 

! * ■ random modulation produces an increase of the cost. We discuss the implications of this phenomenon 

for those overdamped systems whose control parameters are intrinsically noisy, presenting a case 
study based on the example of a Brownian particle optically trapped in an oscillating potential. 
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D 1 I. INTRODUCTION 

S 

Noisy systems are ubiquitous in natural and engineered phenomena. The presence of noise becomes particularly 
■ evident when we move down into molecular-scale phenomena: the thermal noise, responsible for the Brownian diffusion 
' of particles, is omnipresent. However, noise is also intrinsic to many macroscopic systems [l|, [H, [H 0| : stock markets, 

1 t , population dynamics, and traffic flows, all present some degree of noise. 

Even though in the last years the constructive role of noise has been appreciated in many physical, chemical, and 
biological phenomena - examples include stochastic resonance [!, 0, 0, B, 1| , noise induced transitions [13, [H| , noise 
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induced transport [12 . 1 1 31 ] , stochastic resonant damping 14 1 -, there are many situations in which the intrinsic noise 
of a system is still a nuisance that one wants to keep under control, and minimize if possible [ill ] . 

Noise suppression is a crucial task at all scales. Microscopic and nanoscopic phenomena have to deal with thermal 
| noise. Complex pricing systems, such as the Black-Scholes option pricing model [16[, have been developed for dealing 
with the noise present in the stock markets. Given the insolubility of the multi-body problem, noise has to be dealt 
with in the planning of satellites' trajectories. In all these cases, one needs to exert some kind of control on the system 
in order to minimize its intrinsic noise. Often these actions are controlled by some input parameters, which may also 
^/-^ ■ vary over time either deterministically or randomly. 

| Here we will focus on an Ornestein-Uhlenbeck equation, which describes an overdamped system. Such equation 
describes a very wide class of systems, and it has successfully been applied to systems as diverse as macromolecules 



that follow the Hooke's law, Brownian particles, electronic devices, and mesoscopic chemical reactions 17j. As a 
. simple example let us consider the diffusion of a Brownian particle. The particle position variance is reduced if the 
' particle is confined in a potential well. This potential well can be produced by various means: by a molecule that binds 
the particle, by hydrodynamic focusing, or by optical or magnetic tweezers. All these means have in common that 
they exert a restorting force on the particle whenever it is displaced from the desired position. A tighter confinement 
of the particle is achieved by increasing the stiffness of the link, but a higher stiffness implicates undergoing an higher 
• ■ cost to run the system. 

Here we analyze the cost of noise suppression in non-equilibrium systems and introduce a cost function to quantify 
the effort done to control such a system. We show that the minimum cost is achieved when the system control param- 
eters are held constant. We find that any additional deterministic or random modulation of the control parameters 
entails an increase of the cost function. 



II. MODEL 



We consider a dynamical system driven by a Gaussian white random process dB t (Wiener process) , whose state St 
freely evolves according to dst = dB t . By introducing a restoring force, characterized by a constant stiffness k, the 
system can be forced to fluctuate around a state a, with variance <r 2 = l/2k. The stochastic dynamics of the system 
(known as Ornstein-Uhlenbeck process) is described by 

ds t = -k(s t - a)dt + dB t . (1) 
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We can now analyze the effect of fluctuating parameters by letting the mean state at and the stiffness kt - and 
therefore the intrinsic fluctuations <j 2 a = l/2k t - as generic processes independent of dB t and possibly dependent 
between themselves with E[a t ] — and E[k t ] > 0, where E[-] denotes expected value. These two conditions guarantee 
the long term stability of the system. In the case of a Brownian particle, they assure that the particle will not 
eventually escape from the potential. The former condition [E[at] = 0), in particular, signifies that the potential 
keeps on oscillating around the state s = 0; the latter (E[kt] > 0) that the average stiffness is positive. In particular, 
all our conclusions apply to the case in which a t and k t are functions of an Ornstein-Uhlenbeck process. Notice also 
that the conclusions also apply to the case in which at and kt are deterministic function, considering that these are a 
special case of random function. The system time-evolution obeys the equation: 



dst — —k t (st — at)dt + dB t . (2) 

As we have already mentioned, the diffusion of a Brownian particle in a time-varying potential is an example 
of processes described by ((2|) [l4], Gil ■ A free particle diffuses in such a way that the variance of its position grows 
linearly with time. The diffusion process can be partially suppressed by confining the particle in a potential well 
which, for example, can be produced by an optical trap [201 ] . In the presence of an optical trap, whose center and 
stiffness oscillate, the particle position obeys to equation 0, with a t being the center of the trap and k t its stiffness. 

With our analysis we aim to find the output variance of the state of the system described by ([2]) when the parameters 
vary over time in an arbitrary fashion. In particular we will identify four contributions in the total variance: intrinsic 
variance (f^uy stiffness variance <J 2 s ^ k y equilibrium variance <\< a y interplay covariance p s ( a k)- All these cases may be 
studied experimentally. 

Instead of attacking the general case directly, we will proceed by steps, starting by investigating some limiting cases. 
This will permit us to build up the necessary intuition and to gain useful insights into the phenomenon. 

We note that in the cases we study the Ito and Stratonovich approaches to stochastic integration are mathematically 
equivalent, because the diffusion term is constant [21, pp. 35-37]. Here we arc considering the system's steady state, 
but the conclusions apply with little variations also to the transient. 



A. Stationary case (a t = and kt = k) 



The simplest case is when the equilibrium position of the harmonic potential does not oscillate (a t = 0) and its 
stiffness is kept constant (k t = k > 0). This is the benchmark against which all other results will be compared. 
Equation (2) simplifies as 



ds t = —kstdt + dB t . 

Its solution can be find multiplying by the integrating factor e kt and comparing with d ( e^s^j 
The solution is 



(3) 



ke kt s t dt + e kt ds t . 



s t = e kt x +e kt 



l dB„ 
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where the limit has been taken for large t. 

It follows that the mean of the system is E [s t ] — because it is an Ito integral, and its variance is 



E\sj] = E 



-2k(t- 



ds — 



-2k(t-u) du = 



1 - e 



-2kt 



2k 
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where the Ito isometry E 



(Slf{u,u)dB u ] 
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= E 



fof 2 (u,u)du 



has been used [2l|, [2S 



(4) 



(5) 



We can therefore identify the intrinsic variance as a contribution to the total variance of the system: 



1 

2k' 



(G) 



As we will see this is the minimum variance that can be achieved for a given value of the cost function, i.e. for a 
given k. 
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B. Oscillating kt (at = 0) 



When at = 0, equation (2) simplifies as 

ds t = —k t s t dt + dB t , (7) 
where k t is an ltd process independent of B t . Again the solution can be calculated multiplying by the integrating 

r f k u du * . ,,, , / f k u du \ j f k u du ,. . f k u du , T , i , • 

factor e J o and comparing with a e J o s t = k t e J o stdt + e J o dst- Its solution is 



— f k u du , — f k u du f f k v dv 
= P. J a Xq +e Jo / eJa 



st = e 



J dB„ 



J k v dv 



dB„ 



(8) 



where the first term vanishes for large t because E[k t ] = k > 0. 

We can therefore calculate the mean and the variance of the system. 



E [ St ] = E 



J k v dv 



dB„ 



= 0, 



because it is an ltd integral, and 



E[sl\ = I E 



-2 J' k v dv 



du> / e-^-^du = al 



where we have used the ltd isometry and Jensen inequality [23j E 



-2 J* k v dv 



du 



(9) 



(10) 



> e 2£; [ fc *K* w )du integrated over 



time with E[k t ] = k. 

We can now identify the stiffness variance as a contribution to the total system variance, caused by the variation 
of the stiffness: 



7 s{k) 



-2k(t-u) £ 



-2 J t (kv-k)dv ^ 



du. 



(11) 



C. Oscillating at (k t = k) 

The case when the equilibrium position of the potential at is oscillating, while k t = k remains constant, was 
investigated both theoretically and experimentally in [l4| . However, there a different approach was applied and it can 
be useful to obtain the same result expressed in the current formalism. Equation (2) becomes: 



(12) 



dst = —k(st — at)dt + dB t . 

It can again be solved by multiplying by the integrating factor integrating factor e^* and comparing with d ( e kt s t ) = 
ke kt s t dt + e kt ds t . Its solution is 



St = e kt s +ke 



-kt 



e ku a u du + e- kt / e ku dB u >k / e ki - t ' u ^a u du+ / e k{t - u ^dB u 



(13) 



Since the process at is independent of B s , in the calculation of the variance of st the contributions of the two 
integrals can be separated 



s(o)) 



where the equilibrium variance 



< (a) = k z E 



t - 

e k(t ~ u) a„du 



(14) 



(15) 



is the contribution to the variance of the system due to the oscillation of the equilibrium position of the potential. 
The second term is the one corresponding to the stationary state. More details and a discussion of how this effect 
produces the stochastic resonant damping can be found in [14j |. 



4 



D. General case - oscillating a t and kt 

In the general case given by equation (2), again the solution can be calculated multiplying by the integrating factor 



f k u du 
3 Jo 



and comparing with d ( kudu St j = k t e$o k " du stdt + eh kudu dst- The general solution is 



r-t ft ft r-u r-t ft ft r-t ft r-t 

— I k u du , — / k u du I I kvdv, 7 , — / k u du I \ k v dv I I k v dv , , , / / k v dv 

st = e J o so+e J ° / e J o k u a u du + e J o / e J o dB u — > / e J ^ k u a u du+ / e J « 



k u a u du+ I eJu ° dSu- 

For large £ following a procedure similar to the previous cases and finally the variance of the system in the general 
case is given by 



E\s 



E 



^— J k v dv 



dB, 



E 



(J e$u kl ' dv k u a u du^J 



a s(i) + a s(k) + E 



k v dv 



- a2 s(i) + a l(k) + Ps(ak) + t7 s(o ), 



where 



Ps(ak) 



= E 



J k v dv 



a u du 



k 2 E 



t . I 

e k{t ~ u) a u du 



(17) 

(18) 
(19) 

(20) 



is the interplay covariance, which can be either positive or negative. However, the total variance is always larger than 
the intrinsic variance, since, as it can be seen from equation (|19|) . the overall contribution due to the oscillation of the 
stiffness and the equilibrium position is always positive: 



a l(k) + Ps(ak) + 0"s(a) > °- 



(21) 



E. Cost function 



As Eq. pT|) shows that one can use different protocols in order to change the output variance of a given intrinsically 
noisy system by means of the modulation of the system parameters. Now we do the next step and we ask the key 
question of this study: how one can compare the protocols from the point of view an effort applied to change the 
output variance? To deal with such a question mathematically, we suggest to introduce a cost function. 

The idea of a cost function, sometimes referred to as objective function, is very well established in the fields of 
economic optimization [24[ and in engineering (25[: it permits one to compare the performance of systems that work 
under different conditions. Typically for a given cost one looks for the parameters that provide the best performance 
(the smallest variance in our case). We introduce here the idea and the importance of a cost function in the study 
the confinement of overdamped systems. An appropriate cost function needs to describe the overall effort spent in a 
system to achieve its confinement. 

For a stationary system, the stiffness k fully describes the confinement effort. Indeed, as we have seen, the output 
variance of a stationary system is inversely proportional to the stiffness. Therefore to define the cost function as 

C = k seems rather natural. 

We introduce a similar cost function for systems whose parameters arc modulated over time. First, as it was shown 
in [l4| , the modulation of the mean state at does not affect the effort made to confine the system; we therefore need 
to consider only the modulation of the stiffness in order to introduce the cost function. For the systems where the 
stiffness k t varies over time, we suggest to use as a cost function the average value of the stiffness : 

C = E[k t \. (22) 

As seen the cost function of a stationary system calculated using this formula has the same value as it was defined 
before. 

Let us compare the variance of a stationary system and the same system but with modulated parameters, assuming 
that the cost functions are equal for both systems. To maintain the cost function of the system to be constant, we 
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must keep the average stiffness invariant E[kt]=k. From this condition it follows that the intrinsic variance is 
also constant, and, as a straightforward consequence of equation l|19p. it coincides to the minimum variance. This 
means that for a given value of cost function the output variance of the system with modulated parameters is bounded 
by its intrinsic variance: 

E[st] = a%) + cr 2 s(k) + p s(ak) + a 2 {a) > of (i) . (23) 

This can equivalently be stated as the fact that, for a given cost function, any additional deterministic or random 
modulation implicates a larger system variance. For a given value of the cost function the minimum of the variance is 
achieved when the control parameters are constant. From another point of view it means also that for a given system 
variance, any additional deterministic or random modulation produces an increase of the cost function. 

III. BROWNIAN PARTICLE EXAMPLE 

Experimental results related with our study were presented in [l8| where the dynamics of a Brownian particle held 
in an optical trap with modulated position and stiffness was measured. For a given experimental configuration the 
stiffness of the (stationary) trap, and, therefore, the achieved confinement is proportional to the optical power P used 
to create the trap (k oc P). Therefore a higher confinement requires an higher optical power and the cost function of 
the system is defined as C — k cx P. 

In the presence of an optical trap, whose center and stiffness oscillate, the particle position obeys to equation @, 
with a t being the center of the trap and k t — 1/a 2 its stiffness. In this case, the cost function is C = E[k t ] oc E[P(k t )], 
where the optical power P(k t ) that must be used to create the optical trap also fluctuates. 

By using of a specific protocol of modulation of the trap parameters a reduced variance of the observed particle 
position as compared to a stationary trap was observed. We analyzed the experimental data calculating the cost 
function for the stationary and modulated traps. When a stationary trap was used (stiffness k^ = 3.7 pN/fim) the 
output variance of the particle position was a 2 ut — 1087 nm 2 (Fig. 1(b) of [Hj]). With the oscillating trap, the output 
position variance was indeed reduced to <j 2 ut = 764 nm 2 (Fig. 1(c) of Ref. [Hj]). From the data presented in [18j 
we calculated the average stiffness as k = E[k t ] — kpk(k)dk, therefore substituting the time-average with the 
average over pfcffc) the probability density of k t (Fig. 1(b) of Ref. [18]), which can be computed from equations (2) 
and (3) of Ref. The average stiffness of the modulated trap (and therefore the average optical power introduced 
in to the system and the cost function of the process of the noise reduction) was found considerably bigger that in 
the stationary case (k = 6.8 pN/ fim). Therefore, we can conclude that the higher confinement of the particle position 
in the trap was achieved not due to the addition of noise to the trapping parameters, but to the higher average 
trapping power, while the added noise slightly increases the output variance. We notice that a stationary trap with 
such average power and the same experimental configuration would produce even smaller output variance than the 
modulated trap did. 

IV. CONCLUSIONS 

We have shown that to suppress the intrinsic noise of a system entails a cost. For a noisy system in the overdamped 
regime controlled by a fluctuating input parameter, we have found out that the minimum cost is achieved when the 
system control parameters are held constant: any additional deterministic or random modulation of the parameters 
yields an increase of the cost function. It is not possible to reduce the output noise of a system below a threshold 
value, corresponding to constant input parameters without increasing the cost function of the system. This has 
important implications both from the fundamental point of view, in order to understand many natural phenomena - 
for example, the natural optimization of cellular molecular phenomena -, and from the engineering one, where it can 
give a guidance in the management of the intrinsic noise of a system. 
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